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ABSTRACT 

Aims. Reliable timing calibration is essential for the accurate comparison of XMM-Newton light curves with those from other obser- 
vatories, to ultimately use them to derive precise physical quantities. The XMM-Newton timing calibration is based on pulsar analysis. 
However, as pulsars show both timing noise and glitches, it is essential to monitor these calibration sources regularly. To this end, 
the XMM-Newton observatory performs observations twice a year of the Crab pulsar to monitor the absolute timing accuracy of the 
EPIC-pn camera in the fast Timing and Burst modes. We present the results of this monitoring campaign, comparing XMM-Newton 
data from the Crab pulsar (PSR B053 1 +2 1 ) with radio measurements. In addition, we use five pulsars (PSR J0537-69, PSR B0540-69, 
PSR B0833-45, PSR B 1509-58 and PSR B 1055-52) with periods ranging from 16 ms to 197 ms to verify the relative timing accuracy. 
Methods. We analysed 38 XMM-Newton observations (0.2-12.0 keV) of the Crab taken over the first ten years of the mission and 13 
observations from the five complementary pulsars. All the data were processed with the SAS, the XMM-Newton Scientific Analysis 
Software, version 9.0. Epoch folding techniques coupled with^ 2 tests were used to derive relative timing accuracies. The absolute 
timing accuracy was determined using the Crab data and comparing the time shift between the main X-ray and radio peaks in the 
phase folded light curves. 

Results. The relative timing accuracy of XMM-Newton is found to be better than 10~ 8 . The strongest X-ray pulse peak precedes the 
corresponding radio peak by 306±9 ps, which is in agreement with other high energy observatories such as Chandra, INTEGRAL and 
RXTE. The derived absolute timing accuracy from our analysis is ±48 ps. 

Key words, stars: neutron stars - pulsars: individual: PSR B0531+21, PSR J0537-69, PSR B0540-69, PSR B0833-45, PSR B1509- 
58, PSR B 1055-52 - supernova remnants: Crab - X-rays: stars - instruments: EPIC-pn - data analysis: relative timing, absolute 
timing 

1. Introduction ground calibrated with respect to relative timing, but due to a 

,., . . . , ,. „ . . , limited calibration time budget, the end-to-end system for abso- 

A reliable timing calibration is essential for all timing data anal- lute ^ wafj neyer checked Qn ^ nd The ^ 

yses and the physics derived from those. Irregularities in the for ^ ^ ^ Crab wa§ ted tQ haye afl acc 

spacecraft time correlation the on-board instrument oscillators of Ap/p < 1Q _ 8 before launch Foj . ^ absolute ^ a m _ 

or data handling unit and the ground processing and data anal- kement f AT < lms was given. 

ysis software can lead I to errors in relative and absolute infor- XMM-Newton Hansen et all Eoolh was launched in 

mation pertaining to the timing behaviour or astrophysical ob- „ , , „„„ ' ~ ~~z , — ~i, „ , „ 

. . ° . , , . , , December 1999 with an Anane 5 rocket from French Guay ana. 

lects. The timing of the XMM-Newton observatory is evaluated T . . . . . . . „ , ,„, ,. , , 

• xrww *r , i~ T~\Tf~i , • , It operates six instruments in parallel on its 48 hour highly 

using XMM-Newton s EPIC-pn camera that has been extensively „. \. , ... „ T . , . , co . , 

b r J elliptical orbit: three Wolter type 1 telescopes, with 58 nested 
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Grating Spectrometers (RGS; Iden Herder et all 1200 ll) . In 
addition, a 30 cm Ritchey Chretien optical telescop e, the Optical 
Moni tor, is used for optical observations (OM; iMason et al.L 
1200 ll) . EPIC consists of three cameras: the two EPIC-MOS 
cameras use Metal-Oxide Semiconductor CCDs as X-ray de- 
tectors, while the EPIC-pn camera is equipped with a pn-CCD. 
All three have been especially developed for XMM-Newton 
dPfeffermann et al.L 119991: iMeidinger et al.L I1999L iTurner et all 
120011) . 

In this paper we determine the relative timing accuracy of 
XMM-Newton 's EPIC-pn camera using all available observa- 
tions of the Crab pulsar in combination with other isolated pul- 
sars in order to extend our analysis to a broader variety of 
sources. Preliminary results on the relative timing accuracy of 
XMM-N ewton using the Crab pulsar and the other pulsars can be 
found in Caballero et a 

In this work we use only the Crab pulsar X-ray observations 
to determine the absolute timing accuracy. However, as this is 
done in reference to radio timing, it is limited to the accuracy of 
the radio ephemerides. 

We also see the paper as a summary of "how to perform" rel- 
ative and absolute timing analysis with XMM-Newton and what 
timing accuracy the user can expect for different targets. 

This paper is organised as follows. In Sect. |2]we give a de- 
scription of the targets used for the timing evaluation, followed 
by some technical comments on our data analysis in Sect. [3] The 
relative and absolute timing results are presented in Sect. [4] and 
Sect. |5] A short description of the XMM-Newton 's EPIC-pn cam- 
era is given in Appendix lAl 

2. Observations 

All pulsars used in our analysis are isolated. We concen- 
trated primarily on the Crab pulsar (PSR B0531+21) as ra- 
dio ephemerides are provided monthly by the Jodrell Bank 
Observatory^]. The other pulsars have been chosen to include a 
range of periods and pulse profiles, with which to check the rel- 
ati ve timing. Some of these pulsar observations were reported 
bv lBecker et al.l d2002l) as a summary of first results from XMM- 
Newton. 

Tables [T] and [2] summarise the data used and the results 
obtained from all the Crab observations studied and all the 
other pulsars respectively. Column 1 gives the observation ID 
(OBSID) used for identifying XMM-Newton observations, fol- 
lowed by the satellite revolution ("Rev.") in which the obser- 
vation was done, the data mode, and the filter used. Column 5 
indicates whether the observation is affected by telemetry gaps 
(due to a full science buffer), and column 6 gives information on 
time jumps during the observation (see the footnote of the table 
for explanation). Column 7 lists the start times ("Epoch") of the 
observations in MJD, followed by the exposure ("Obs. Time") 
in ks. Columns 9 and 10 list the pulse periods of the Crab pulsar 
in the radio at the time of the XMM-Newton observations (in- 
terpolated using the information provided by the Jodrell Bank 
Observatory) and the measured X-ray period, respectively. Red. 
X 2 (column 11) gives the reduced^ 2 values found at the maxi- 
mum of the respective x 2 distribution of the period search (the 
number of degrees of freedom, dof, was always 100 for the Crab 
pulsar), and "FWHM" is the full width at half maximum of the 
\ 2 distribution. AP/P is the relative difference between the radio 
and the X-ray period (Eq.Q]in Sect. 1311 The "Phase Shift" (last 
column) shows the measured time shift of the main peak in the 



pulse profile between the X-ray and radio profiles, as explained 
in Sect. 13.21 All uncertainties given are at the 1 <x (68%) level. 
The ephemerides of all the targets used in the analysis are shown 
in Table|4] Fig.Q]shows the pulse profiles for all the pulsars anal- 
ysed in this paper. 

2. 1 . The main XMM-Newton timing monitoring source: PSR 
B0531+21 (The Crab pulsar) 
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Since the discovery of the Crab Pulsar (Staeli n & ReifensteinL 
1968), the Crab has been one of the best studied objects in the 
sky and it remains one of the brightest X-ray sources regularly 
observed. As a standard candle for instrument calibration, the 
33 ms Crab pulsar has been repeatedly studied (monitored) by 
many astronomy missions in almost every energy band of the 
ele ctromagnetic spectrum. Ho wever, recent analysis presented 
bv lWilson-Hodge et alJ (1201 lh showed that the flux of the Crab 
is not constant on long timescales at high energies. These flux 
variations seem to be related to the nebula and correspond to a 
flux drop of ~7 % (70 mCrab) over two years (2008-2010). This 
might affect the status of the Crab as a standard candle in the 
future. 

In the X-ray regime its pulse profile exhibits a double peaked 
structure with a phase separation of 0.4 between the first (main) 
and the second peak. X-ray emission at all phases, incl uding 
the pulse minimum, was discovered bv lTennant et al.l (1200 lh us- 
ing the Chandra observatory. Measurements of X-ray to radio 
delays between the arrival times of the main pulse in each en- 
ergy range of the Crab pulsar have been reported using all high- 
energy instruments aboa rd INTEGRAL dKuiper et all 20031) and 
RXTE dRots et alll2004b . The time delays were determined to be 
280±40 fis and 344+40 //s respectively. 

The Crab pulsar has been observed bianually to monitor the 
timing capabilities of XMM-Newton. Over the years an observa- 
tion strategy has been established that makes very efficient use of 
the limited calibration time budget. XMM-Newton generally ob- 
served the Crab pulsar three times per orbit for 5 ks: at the begin- 
ning, in the middle, and the end of that orbit. These campaigns 
were carried out in spring and autumn when XMM-Newton has 
a different location in its orbit with respect to the Sun-Earth sys- 
tem. This guarantees the monitoring of the dependency of the 
timing with respect to XMM-Newton 's orbital position. Eventual 
irregularities in relative timing with respect to the orbital posi- 
tion could then be identified. A total of 38 observations with ex- 
posure times between 2 ks and 40 ks have been analysed in this 
paper. See Table[T]for details of these observations. 

2.2. Other useful pulsars for relative timing analysis 
2.2.1. PSR J0537-69 

PSR J0537-69 is a young pulsar, about 5000 years in age, lo- 
cated in the Large Magellanic Cloud. It is embedded in the 
supernova remnant N157B and is considered to be the oldest 
known Crab-like pulsar. It is a v ery fast-spinnin g pulsa r with 
a period of 16 ms, discovered by Marshall et al. (1998) using 
RXTE. No significant rad io signal above a 5cr th reshold has been 
detected from the pulsar (ICrawford et al.L l2005h . In the X-ray en- 
ergy range, RXTE has monitored PSR J05 3 7-69 f or seven years 
(Mar shall et all I2004L iMiddledit ch et all 120061). providing a 
compl ete study of the behaviour of the pulsar. Middleditch et al.l 
(2006) reported 23 sudden increases in frequency, called glitches 
and present in most of young pulsars. Due to this highly ir- 
regular activity (a glitch every ~4 months) a contemporaneous 
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ephemeris is important. Its pulse profile in the X-ray regime is 
characterized by a single narrow peak. 

See Table|2]for details of the observations. Our 36 ks obser- 
va tion coincides with the RX TE monitoring campaign presented 
by Middleditch et al. (2006) and a good ephemeris was therefore 
guaranteed. 



2.2.2. PSR B0540-69 

Th is young pulsar (~1 500 years) was discovered in soft X-rays 
bv lSewardetaTl(ll984l) with a period of 5 1 ms, in the field of the 
Large Magellanic Cloud and it is considered to be a Crab-like 
pulsar. Its pulse shape does not appear to chan ge significantly 
from optical to hard X-rays (Ide Plaa et all "2003). The pulsed ra- 
dio emissi on was discovered in late 1989 appearing as a faint 
source (Manc hester et all 1 19931) and presenting a complex pro- 
file, very different from the simple sinusoidal one see n in X-rays 
(FiglD- A glitch was reported by Zhan g et ail (1200 1|) before the 
XMM-Newton observations and confirmed by Living stone et all 
(2005b) using a 7.6 year RXTE campaign. The glitch activity of 
PSR B0540-69 is known to be less than that of the Crab pulsar 
dLivingstone et all l2005bl) but the presence of con siderable tim- 
ing noise was reported bv ICusumano et al.l (120031) using ASCA, 
BeppoSAX and RXTE observations made over a time interval of 
about 8 years. Therefore, despite the low glitch activity, long ex- 
trapolations of the ephemeris might not be reliable. See Table [2] 
for details of the observations. 



2.2.3. PSR B0833-45 (Vela pulsar) 

The Vela pulsar with a period of 89 ms was discovered by 
lLarge. Vaughan & Millsl (f 19681) and it is associated with a 
supernova remnant. It is, with the Crab pulsar, one of the 
most active young/middle-age pulsars known, showing regular 
glitches. These glitches have been intensively studied for the 
Vela pulsar, where a dozen events in different energy ranges 
have been recorded and analysed over th e past three decades 
dHelfand et all I200U [Dodson et all 120071) . Due to these impor- 
tant irregularities close radio ephemerides are needed. 

No Vela timing mode observations have been performed with 
XMM-Newton, but since the period is 89 ms the data in Small 
Window mode (time resolution of 5.7 ms) can be used for our 
purposes. Thus we have analysed the four observations listed in 
Tabled 



2.2.4. PSR B1 509-58 

This young pulsar (~1700 years) is one of the most energetic 
pulsars known and has a pulse period of ~ 15 1 ms. It is associ- 
ated with the supernova remnant G320.4-1.2 and it has been well 
studied in all wavelength s since it was disco vered in the soft X- 
ray band using Einstein (Se ward et all 1 1 98*21) . The pulse profile 
in X-rays of appears to be much broader than in radio, changing 
from a narrow peak shape into a more sinusoidal shape at high 
energies. Monitored by RXTE since its launch, and covering a 
21 year time interval and in conjunction with radio data from 
the MOST and Pa rkes observatories , a detai led timing study has 
been carried out (Livin gstone et all l2005al) . but no glitch was 
found in the entire data sample. This result makes PSR B1509- 
58 probably the only known young pulsar that does not present 
any glitches over long periods of time. This property means 
that it is well adapted to extrapolation over long time intervals 



and useful for absolute timing analyses dRots et all fl998). See 
Table|2]for details of the observations. 

2.2.5. PSR B1 055-52 

PSR B 1055-52, one of the Three Musketeers together with PSR 
B0656+14 and Geminga, is a middl e-aged pulsar with a pe- 
riod of 197 ms. It was discovered by Vau ghan & Largel (1 19721) 
but it was only in 19 8 3 that X-ray emission was first detected 
bv ICheng & Helfandl ([19811) using the Einstein Observatory. 
lOgelman & Finlevl d 19931) detect ed sinusoidal pulsat ions in X- 
rays up to 2.4 keV. More recently, iDe Luca et al.l (120051) showed 
using XMM-Newton data that the pulsed emission is detectable 
up to 6 keV. Most middle-aged pulsars like PSR B 1055-52 
show reduced timing noise and fewer glitches compared to 
younger ones. See Table [2] for details of the observations. 
Results concerning these data have been originally published by 
IDe Luca et all d2005l) . 
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Fig. 1. XMM-Newton pulse profiles of the different pulsars anal- 
ysed. From top to bottom: PSR J0537-69 (obs. ID: 01 13020201), 
the Crab pulsar (obs. ID: 0122330801), PSR B0540-69 (obs. 
ID: 0125120201), the Vela pulsar (obs. ID: 0111080201), 
PSR B 1509-58 (obs. ID: 0312590101) and PSR B 1055-52 (obs. 
ID: 0113050201) with periods of 16 ms, 33 ms, 51 ms, 89 ms, 
151 ms and 197 ms, respectively. The energy band in all cases is 
0.2-12 keV. 
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Table 1. Individual observations from the Crab monitoring. (A description of the columns is given in the text.) 



OBSID 


Rev. 


Mode 


Filter 


Counting mode 


Time lump 


Epoch 


Obs. Time 


Radio period 


X-ray period 


Red. x 1 


FWHM 


AP/P 


Phase Shift 














[MID] 


[ks] 


[ms] 


[ms] 




10"S 


10" 8 


[fXS] 


0122330801 


56 


TI 


T 


NO 


S 


51632.83 


22.58 


33.5083819072 


33.5083817(1) 


433 


1.37 


-0.6±0.4 


-148±2 


0135730701 


234 


B 


T 


NO 


N 


51988.64 


10.00 


33.5213091402 


33.5213091(5) 


197 


3.24 


0.02+0.9 


-337±4 


0153750201 


411 


B 


M 


YES 


N 


52340.68 


4.63 


33.5341004722 


33.5341001(7) 


79 


6.78 


-0.9±2 


-287±4 


0153750301 


411 


B 


M 


YES 


N 


52341.29 


9.00 


33.5341224876 


33.5341223(9) 


143 


3.52 


-0.3+1 


-296+4 


0153750401 


411 


TI 


T 


YES 


N 


52341.42 


9.00 


33.5341272691 


33.5341278(9) 


269 


3.57 


2±1 


-217±2 


0153750501 


411 


B 


M 


YES 


N 


52341.84 


9.00 


33.5341426021 


33.5341426(5) 


136 


3.55 


0.2±1 


-257±4 


0160960201 


698 


TI 


T 


YES 


N 


52913.51 


28.26 


33.5549129334 


33.5549128(4) 


912 


1.58 


-0.3±0.5 


-174±2 


0160960301 


700 


TI 


T 


YES 


N 


52918.30 


10.17 


33.5550871157 


33.5550881(2) 


313 


4.58 


3±1 


-292±3 


0160960401 


874 


B 


T 


NO 


N 


53264.10 


14.72 


33.5676516600 


33.5676516(5) 


250 


2.20 


-0.04±0.7 


603±3 


0160960601 


874 


B 


T 


NO 


N 


53265.57 


3.98 


33.5677052394 


33.5677050(1) 


68 


8.01 


-0.7+2 


8839+5 


0160960701 


955 


B 


T 


YES 


N 


53425.72 


8.20 


33.5735256348 


33.5735259(2) 


133 


3.86 


0.9±1 


-334±3 


0160960801 


955 


B 


T 


YES 


N 


53426.27 


5.00 


33.5735454128 


33.5735452(4) 


80 


6.52 


-0.5±2 


-336±5 


0160960901 


955 


B 


T 


YES 


N 


53426.77 


8.20 


33.5735638700 


33.5735642(0) 


127 


3.90 


1±1 


-357±4 


0160961001 


1048 


TI 


T 


YES 


N 


53612.63 


4.99 


33.5803179048 


33.580317(7) 


178 


11.18 


-0.6±3 


-257±3 


0160961101 


1049 


B 


T 


YES 


N 


53613.14 


5.00 


33.5803366674 


33.5803368(8) 


59 


6.60 


0.6±2 


-335±4 


0160961201 


1049 


TI 


T 


YES 


N 


53613.92 


4.99 


33.5803650194 


33.580364(2) 


164 


10.25 


-3+3 


-286±2 


0160961301 


1140 


B 


T 


NO 


N 


53794.72 


5.00 


33.5869348775 


33.5869347(5) 


90 


6.77 


-0.4±2 


-284±4 


0160961401 


1140 


TI 


T 


YES 


N 


53795.44 


4.99 


33.5869607624 


33.5869607(3) 


262 


6.18 


-0.1 ±2 


-209±2 


0160961501 


1140 


B 


T 


NO 


N 


53796.09 


6.70 


33.5869844395 


33.5869844(8) 


129 


4.47 


0.1+1 


-305+4 


0312790101 


1138 


B 


T 


YES 


E 


53791.22 


42.50 


33.5868077617 


33.5868078(5) 


422 


1.25 


0.3±0.4 


-268±2 


0312790201 


1138 


B 


T 


NO 


N 


53791.75 


5.03 


33.5868269683 


33.5868267(9) 


72 


6.46 


-0.5+2 


-349±5 


0312790401 


1138 


B 


T 


NO 


N 


53791.95 


7.80 


33.5868339621 


33.5868344(5) 


52 


4.19 


1 + 1 


-305+6 


0412590101 


1249 


B 


T 


NO 


N 


54012.06 


6.40 


33.5948315461 


33.5948316(2) 


1 1 1 


4.82 


0.2±1 


-369+4 


0412590201 


1249 


TI 


T 


YES 


N 


54012.72 


5.04 


33.5948556935 


33.5948540(4) 


282 


9.77 


-5+3 


-284±2 


0412590301 


1249 


B 


T 


NO 


N 


54013.38 


8.85 


33.5948794078 


33.5948792(9) 


148 


3.39 


-0.4±1 


-389+3 


0412590601 


1325 


TI 


T 


YES 


N 


54164.32 


5.04 


33.6003636968 


33.6003640(7) 


241 


6.14 


1±2 


-524±2 


0412590701 


1325 


B 


T 


YES 


N 


54164.98 


8.85 


33.6003876860 


33.6003877(0) 


156 


3.50 


0.03+1 


-601+3 


0412591001 


1414 


B 


T 


NO 


N 


54341.11 


13.20 


33.6067862758 


33.6067861(8) 


220 


2.33 


-0.3±0.7 


-336±3 


0412591101 


1414 


TI 


T 


YES 


N 


54341.96 


5.00 


33.6068170822 


33.6068166(3) 


308 


8.62 


-1±3 


-303±3 


0412591201 


1414 


B 


T 


NO 


S 


54342.40 


13.41 


33.6068330929 


33.6068330(5) 


228 


2.40 


-0.1+0.7 


-360±3 


0412591401 


1504 


B 


T 


YES 


N 


54520.72 


21.70 


33.6133108411 


33.6133108(9) 


394 


1.45 


0.1+4 


-357±2 


0412591501 


1504 


TI 


T 


YES 


N 


54521.46 


7.64 


33.6133378231 


33.6133377(7) 


521 


4.53 


-0.1±1 


-282±2 


0412591601 


1504 


B 


T 


YES 


N 


54521.85 


20.20 


33.6133519405 


33.6133519(9) 


382 


1.51 


0.1 ±0.4 


-363±2 


0412591901 


1600 


B 


T 


YES 


N 


54712.22 


5.80 


33.6202664839 


33.6202667(5) 


133 


5.50 


0.8±2 


-282±3 


0412592001 


1600 


TI 


T 


YES 


N 


54712.78 


6.04 


33.6202865236 


33.6202880(2) 


544 


7.19 


4±2 


-207±2 


0412592101 


1600 


B 


T 


YES 


N 


54713.35 


15.39 


33.6203076395 


33.6203076(3) 


365 


1.98 


-0.02±0.6 


-263±2 


0412592401 


1687 


B 


T 


YES 


N 


54885.70 


18.90 


33.6265670615 


33.6265671(0) 


355 


1.67 


0.1±0.5 


-379±2 


0412592501 


1687 


TI 


T 


YES 


N 


54886.52 


4.54 


33.6265967322 


33.6265966(8) 


341 


6.78 


-0.1+2 


-305+2 



Notes. N: no time jump, S: Time jump corrected by the SAS, E: Time jump not corrected by the SAS, data had to be partly excluded from the analysis. TI: Timing Mode, B: Burst Mode, M: 
Medium Filter, T: Thick Filter. Epoch: the epoch of the 1 st XMM-Newton event. Radio period: the radio period extrapolated to the XMM-Newton epoch. The accuracy of the radio periods can be 
considered to be good to 10~ 13 , which is almost negligible compared to the X-ray errors of the period. 



Table 2. Individual observations of the other objects. 



Object 


OBSID 


Rev. 


Mode 


Filter 


Epoch 


Obs. Time 


Reference period 


Reference Epoch 


X-ray period 


Red.* 2 


FWHM 


AP/P 












[MJD] 


Da] 


[ms] 


[MJD] 


[ms] 




10~ 8 .s 


10~ 8 


pep ins^l x 


01 1 1070701 


357 


Tl 


M 


S7717 OS 


JJ.7J 


16 17798071% 


S7760 40796 


16 1779807H^ 


24 


18 


_n oi + 7 


PSR R0540* 


0175170701 


85 


Tl 


M 


51691 54 


17.05 


50 5197855014 


51686 70757 


50 519784C6'! 


99 


12.60 


-2.40 + 12 


PSR R0540* 


0175170701 


85 


Tl 


M 


51691 54 


17.05 


50 5 1 9443 1419 

Jv.J 1 7 II J It^I 7 


57857 86600 


50 519784C6'l 


99 


12.60 


-31.44 + 12 


PSR R0540* 


041 31 80701 


1248 


Tl 


M 


54010 06 


12.64 


50 61 6078661 7 

JVJ. VJ 1 UUZ.OUU 1 z. 


57857 86600 


50 61 5775f9~l 

JU.U1 JZ^Z.J^7^ 


112 


17.34 


-1580 + 17 


rorv Dujtu 


041 ^1 80^01 


1248 


Tl 


M 




15.24 


SO fSl 605Q40R^ 


^78^7 86600 


so (Si s?s?n , i 

jU.UlJiJi^/ 7 


1 16 


14.75 


_ i 504 + i <; 
1J74 in i j 


Vela* 


01 1 1080101 


180 


SW 


M 


5 1 880 00 

J 1 OOVJ.VJU 


37.90 


89 3318630836 

07.JJ lOUJUOJU 


57408 00000 

JZ.TUO .\J\J\J\J\J 


07.JJ lOlO^J ) 


18 


7.90 


-17.35 + 3 


Vela* 


0111080201 


180 


SW 


M 


51880.51 


58.60 


89.3318685573 


52408.00000 


89.331853(2) 


25 


6.01 


-17.84 ±2 


Vela* 


0111080301 


180 


SW 


M 


51881.23 


1.94 


89.3318762404 


52408.00000 


89.33181(3) 


2 


81.00 


-79.22 ± 30 


Vela* 


0153951401 


1169 


SW 


M 


53852.54 


120.87 


89.3531575271 


53193.00000 


89.352786(3) 


17 


8.47 


-415.46 ± 3 


PSR B 1509* 


0128120401 


137 


Tl 


M 


51794.26 


9.60 


151.1145869352 


50352.00000 


151.11357(1) 


198 


190.19 


-673.21 ±25 


PSR B 1509* 


0312590101 


1136 


Tl 


M 


53786.98 


31.95 


151.3773132754 


53385.00000 


151.37723(1) 


694 


58.35 


-52.25 ± 8 


PSR B 1055* 


0113050101 


186 


Tl 


M 


51892.93 


19.90 


197.1118094698 


50256.00000 


197.1118(1) 


30 


151.09 


12.83 ±55 


PSR B 1055* 


0113050201 


187 


Tl 


M 


51893.73 


52.70 


197.1118098756 


50256.00000 


197.11183(5) 


34 


65.81 


11.09 ±24 



Notes. M: Medium Filter; TI: Timing Mode; SW: Small Window Mode. '*': reference ephemeris from radio data; m : reference ephemeris from RXTE data. 
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3. Data analysis 

The data sets were processed using the XMM-New ton Scientific 
Analysis Software, SAS 9.0 (iGabriel et all 12004 . Event times 
were corrected to the solar system bary centre using the SAS tool 
barycen. 

3.1. Relative timing data analysis 

We define the relative timing accuracy as the difference between 
the period measured with XMM-Newton and the period mea- 
sured at radio wavelengths evaluated at the epoch of the X-ray 
observations. This difference is normalised to the pulse period 
measured in radio. 

„ , . . Px-ray(Px-ray) _ Pradio(Px-ray) AP 

Rel. timing := = — (1) 

PradioUx-ray) P 

where 

Px-ray '■ period derived from XMM-Newton 
Pradio ■ period extrapolated from radio ephemeris 
Px-ray : time of the first X-ray event of the XMM-Newton 
observation [MJD] 

We determined the period of the Crab pulsar in X-rays using 
the epoch folding software XRONOS^. The closest available radio 
ephemeris (supplied by the Jodrell Bank Crab Pulsar Monthly 
Ephemeris) before and after the X-ray observation were used to 
interpolate the radio period P for the time of the first X-ray event 
of the XMM-Newton observation in MJD. The interpolated radio 
periods are then used as an initial trial value for the epoch fold- 
ing. The period derivative P provided by Jodrell Bank is taken 
into account when doing the folding of the X-ray data. All rele- 
vant initial and final values are listed in Tables [TJ|2][3]and|4] All 
X-ray pulse profiles shown in Fig. [TJhave been produced using 
the best fit X-ray period. 

The detailed steps of our data reduction are presented below 
in order to provide an example for proper XMM-Newton relative 
timing data analysis. 

1. calibrate the XMM-Newton event list using the SAS routine 
epprocQ 

2. perform barycentre correction using precise coordinate s with 
the SAS routine barycerQ 

3. extract soured 

4. extrapolate the radio ephemeris 

5. period search using ef search from XRONOS (see Table [3) 
which gives the x 1 against the period 

6. period determined through a weighted mean of all values 
within 65% of the ef search^ 2 maximum 

The number of phase bins per period (nphase) in each pulse 
profile was chosen such that the count rate uncertainties in each 

2 XRONOS is part of the HEARSAC software 
(http : //heasarc . gsf c . nasa . gov). 

3 command line set up of epproc: timing=YES burst=YES sr- 
cra=83.633216667 srcdec=22.014463889 withsrccoords=yes 

4 command line set up of barycen: withtable=yes ta- 
ble='bary.ds:EVENTS' timecolumn='TIME' withsrccoordinates=yes 
srcra='83.633216667' srcdec='22.014463889' processgtis=yes 
time=0 

5 Detector coordinates used in the extraction process a) timing 
mode: (RAWX,RAWY) IN box(35,101, 12,100,0), b) burst mode: 
(RAWX,RAWY) IN box(35., 71.5,20,70,0) 



Table 3. Settings for the epoch folding using ef search. 



Object 


nper 


nphase 


dres 


P 


PSRJ0537 


5000 


15 


io- 10 


5.1815 10~ 14 


Crab 


5000 


100 


10 -io 


0.4205 10~ 12 


PSR B0540 


10000 


20 


10 -io 


4.78907 10- 13 


Vela 


2000 


30 


io- 9 


1.25008 10" 13 


PSR B 1509 


50000 


50 


10 -io 


1.53085 10~ 12 


PSR B 1055 


40000 


14 


io- 10 


5.8354873 10~ 15 



Notes, nper: number of periods over which the search is carried out; 
nphase: resolution of the trial folded light curves (number of bins in 
each period); dres: resolution for the period search in seconds; P: period 
derivative. 



bin (determined using the Poisson error on the count rate per 
bin normalised by the bin size) are, on average, not bigger than 
approximately 10% of the total count rate variation in the pulse 
profile of the shortest observation for each pulsar. This value, de- 
termined for each pulsar, is used for all the observations of that 
object. In this way the signal to noise in each bin is sufficient 
to reliably determine any 'smearing out' of the pulse profile due 
to the use of an inaccurate period/ephemeris, essential for deter- 
mining the relative timing precision, as described in Sect. [4] 

3.2. Absolute timing data analysis 

The XMM-Newton EPIC-pn absolute timing accuracy was de- 
termined using only observations of the Crab. The ephemeris 
(epoch, P, P, P) of the nearest radio observation from the Jodrell 
Bank Observatory was used as a reference to obtain the phase 
shift between the time of arrival of the main peak in the X-ray 
profile and the time of arrival of the main peak in the radio pro- 
file, as described in Eq.f2] The phase shift was then multiplied by 
the corresponding X-ray period found during the relative timing 
analysis, as shown in TableQ] 

Phase Shift [// s] := T 0x . ray - T 0radio (2) 

where 

Pox-ray : Time of arrival of the main peak of the X-ray profile 
Poradio : Time of arrival of the main peak of the radio profile 

The phase of the main X-ray peak was determined using a 
pulse profile with 1000 phase bins which was then fitted with 
an asymmetrical Moffat function. The explicit formula for the 
Moffat function is given in Appendix C. We also demonstrate 
how its shape varies when different parameters are modified. 
Fig [2] shows an example of how the phase of one Crab X-ray 
pulse profile (obs. ID: 0122330801) is slightly shifted in phase 
with respect to the radio phase (shown as a red line). 

The following steps describe an example of the data reduc- 
tion carried out on the XMM-Newton data in order to assess the 
absolute timing precision: 

Steps 1-3 are the same as described in Sect. 13.11 

4 fold the X-ray data on the radio period 

5 fit the X-ray pulse profile with a Moffat function 

6 determine the shift between the radio phase zero and the X- 
ray peak. 
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Table 4. Radio and RXTE ephemerides used in the analysis. 



Object 


RA 


Dec 


epoch E 


y(E ) 


HE ) 


KEo) 




(J2000) 


(J2000) 


[MJD] 


[mHz] 


IxlO-' Hz^ 1 ] 


[xlO- 21 Hz.?" 2 ] 


PSR J0537 1 


05:37:47.20 


-69:10:23.00 


52260.40296 


62.02279999(0) 


-1.993398(1) 


6.80(0) 


Crab 2 


05:34:31.97 


22:00:52.07 










PSR B0540 3 


05:40:11.22 


-69:19:54.98 


51686.20757 


19.794507(5) 


-1.881021(0) 


3.79(9) 


PSR B0540 4 


05:40:11.22 


-69:19:54.98 


52857.86600 


19.77552961(8) 


-1.87288(3) 


4.30(2) 


PSR B0540 4 


05:40:11.22 


-69:19:54.98 


52857.86600 


19.775529611(3) 


-1.872853(1) 


4.18(1) 


Vela 5 


08:35:20.61 


-45:10:34.87 


52408.00000 


11.193503640388(1) 


-0.156027(3) 


0.64(1) 


Vela 6 


08:35:20.61 


-45:10:34.87 


53193.00000 


11.192447207118304(3) 


-0.155502(8) 


0.52(7) 


PSRB1509 7 


15:13:55.62 


-59:08:09.00 


50352.00000 


6.625918674074(4) 


-0.673579(0) 


1.95(0) 


PSRB1509 8 


15:13:55.62 


-59:08:09.00 


53385.00000 


6.608333867737907(3) 


-0.66852(2) 


1.91(4) 


PSRB1055 9 


10:57:58.84 


-52:26:56.30 


50256.00000 


5.073283989285(9) 


-0.015019(5) 


3.24(0) 



Refer ence s. (1) dMiddleditch et all 120061): (2) Jodrell Ba nk C rab Pulsar Monthly Ephem eris; (3) dCusumano et all l2003h : (4) dJohnston etafl. 
120041) : (5) dRomani et all 120051) : (6) dDodson et all 120071) : (7) dLivingstone et all 120053) : (8) Parkes Observatory (private communication); (9) 
Princeton Database. 



3500 7 
3000 : 



2500 r 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 
Phase (set at radio epoch0=0) 



Fig. 2. Crab pulse profile (obs. ID: 0122330801) with 1000 
phase bins which was used to determine the phase of the main 
peak. The phase of the main peak at radio wavelengths is shown 
as a solid red vertical line. 



3.3. Evaluating the efficiency of automatic corrections made 
to event time jumps by the SAS 

In order to do proper timing analysis with the EPIC data, ev- 
ery event detected on board has to be assigned a correct pho- 
ton arrival time. The transformation from readout sequences 
by the EPIC camera to photon arrival times of each photon 
is performed by th e EPEA (European Photon Event Analyser, 
iKuster et all fl9 99). The absolute timing adjustment from On 
Board Ti me to UTC is don e with the XMM-Newton time cor- 
relation (Kirsch et al., 2004). A hardware problem in the EPEA 
can produce time jumps in some observations, which have to be 
corrected. 

A time jump can affect the timing accuracy by broadening 
the ^-distribution d uring the epoch fo lding search or by pro- 
ducing 'ghost peaks' dKirsch et alll2004h . 

To find time jumps reliably one needs to have as accu- 
rate CCD frame times as possible. The fram e times for all 
EPIC-pn modes using all available archive data (Frevber g et all 
l200l were remeasured and cross-checked with an independent 
method, analysing the evolving differences between consecutive 



events. If these differences are correct, the slope of the relation 
between the difference between the arrival times of consecutive 
events and the same quantity modulo the frame time is zero. 
Fig. [3] shows the effect of the frame time recalibration, which 
brought the values of this slope very close to zero for all ob- 
servational modes. Constant 'timediff' values' indicate constant 
frame time. 

Recalibration of the time jump detection algorithm of the 
XMM-Newton SAS has been done with the refined frame times. 
A search for time jumps in all available XMM-Newton archive 
data up to revolution 1061 showed a significant number of time 
jumps in the data for each observational mode. The applica- 
tion of the new algorithm reduced the number of remaining time 
jumps significantly. The effects of the SAS JUMP .TOLERANCE pa- 
rameter in the new algorithm (Versions 2-7) are shown in Fig. [4] 
for each EPIC-pn mode. The Version shown in Fig. |4]repre- 
sents the remaining time jumps with the old frame times and the 
old algorithm and it is included for comparison. Version 1 was 
obtained using the new frame times but with the old algorithm. 
While the rate of non-corrected time jumps (averaged over all 
pn-modes) was 2.8 per 100 ks before the implementation of the 
SASv8.0 time jump correction, just 0.3 time jumps per 100 ks 
remained uncorrected after its implementation. A breakdown of 
time jumps for each EPIC-pn mode is given in Table [5] This 
new improved algorithm has been implemented in the SAS as 
the default setting since version 8.0 (SAS_JUMP_TOLERANCE = 
22.0). Table[T]indicates for which Crab observation a time jump 
has been corrected, and where data had to be excluded from the 
analysis. 

In order to identify possible remaining time jumps, the 
data can be processed without the "fine-time"-correction, 
i.e., epframes set="infile_pn" eventset=events.dat 
gtiset=tmp_g . dat withfinetime=N. Then the time At 
between successive events is calculated and divided by the 
frame time, FT, of the relevant mode (FF Mode: 73.36496 ms, 
eFF Mode: 199.19408 ms, LW mode: 47.66384 ms, SW Mode: 
5.67184 ms, Timing Mode: 5.96464 ms, Burst Mode: 4.34448 
ms). 

A time jump is shown to exist when At/FT is different from 
an integer by a quantity larger than a tolerance parameter. Only 
those time jumps which happen to be an integer multiple of the 
relevant FT would not be found with this method. It is important 
to notice that the tolerance acceptable between At and the full 
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Fig. 3. The time difference between consecutive events deter- 
mined using the slope of a linear fit to the time differences 
of consecutive events modulo the frame time. This is plotted 
against time. Upper panel with old frame times, lower panel 
with the refined ones. The different modes of EPIC-pn are rep- 
resented by different colours: Full Frame (black, Extended Full 
Frame (dark blue), Large Window (green), Small Window (red), 
Timing (light blue), Burst (yellow). 



frame time should not be bigger than (20/48828. 125 xFT). 



Table 5. Mean rate of residual uncorrected time jumps per 
ks. 



100 



pn instrument mode 


Pre SASv8.0 


Post SASv8.0 




per 100 ks 


per 100 ks 


Full Frame (FF) 


0.6 


0.2 


Extended Full Frame (eFF) 


2.1 


0.2 


Small Window (SW) 


3.7 


0.7 


Large Window (LW) 


3.9 


0.1 


Timing 


1.2 


0.1 


Burst 


5.3 


0.3 


Overall mean 


2.8 


0.3 



$ 30 



E 

3 



I 20 



C 




Fig. 4. Remaining time jumps in the EPIC-pn data after all 
SAS correction algorithms. The triangles correspond to the re- 
maining time jumps (in percent) for each Epic-pn mode: FF 
(black), eFF (blue), SW (red), LW green), Timing (light blue) 
and Burst yellow). Filled circles represent the overall remaining 
time jumps for all XMM-Newton observations up to revolution 
1061. The numbers stand for the different processing versions 
of the algorithm. Versions: 0: old frame times (oft) and old al- 
gorithm, 1: new frame times (nft) and old algorithm, 2: nft and 
SAS_JUMP_TOLERANCE (SJT) 19.0, 3: nft and SJT 20.0, 4: nft and 
SJT 21.0, 5: nft and SJT 22.0, 6: nft and SJT 23.0 and 7: nft and 
SJT 24.0 



4. Relative Timing accuracy of XMM-Newton 

The relative timing accuracy of the EPIC-pn camera has been 
studied using all six pulsars (see Fig. |TJ», presented in Sect. [2] 

As described in Appendix B, the FWHM of the y 2 curve 
obtained during the period search analysis can be expressed 
in terms of the period and the exposure time of the observa- 
tion (Eq. IB.2l ). From t his expression, and usin g the Independent 
Fourier Space (IFS; Ide Ofia Wilhelmil 120031) . approach dis- 
cussed in the appendix, an empirical formula for the error on 
the X-ray period was found (Eq. [3J. 



5P = 



FWHM 

~doJ 



(3) 



where dof is the number of degrees of freedom (number of phase 
bins used to construct the pulse profiles minus the number of 
variables). The number of bins used in the pulse profiles are 
shown in Table [3] 

The relative timing accuracy was defined by Eq. [T]and there- 
fore its error will depend mostly on how accurately the radio 
and X-ray periods can be measured. Other factors that could af- 
fect the relative timing accuracy are discussed in Appendix B. 
Considering that the radio period measurements are more accu- 
rate than the X-ray periods (usually by 1-3 orders of magnitude; 
however, DM variations can sometimes cause problems and the 
time resolution of the radio telescopes has to be monitored), it 
was assumed in our analysis that their errors were negligible 
compared to the error on the X-ray period. Thus, it can be found 
that the relative timing, AP depends exclusively on the error of 
the X-ray period, 6P as shown in Eq. [4] 



AP * 6P 



(4) 
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The relationship described in Eq. |4] allows a goodness of 
fit study of our measured AP compared to the empirical 6P de- 
scribed in Eq. [3] The empirical 6P was considered as an upper 
limit on an accurate AP measurement. A comparison between 
the observed relative timing accuracy in absolute value and nor- 
malised by the corresponding period, |AP/P| (symbols) and its 
"expected" value obtained from Eq. [3] and Eq. [4] (lines) is pre- 
sented in Fig. |5] as a function of the exposure time and in Fig. [6] 
as a function of date. 



1(TF 




10 2 1 3 1 4 1 5 

Time Base [s] 



Fig. 5. The relative timing accuracy of XMM-Newton 's EPIC-pn 
camera. Comparison of the expected |AP/P| using the assumption 
that the relative timing depends exclusively on the error of the X- 
ray period. 6P can be empirically expressed as shown in Eq. [3] 
and |4] (given as the lines in the Figure) and the measured one 
obtained by comparing our X-ray period with the extra-inter- 
polated radio period (symbols). 

For observations where we have very reliable radio /X-ray 
ephemerides, the observed data points are below the lines of the 
estimated accuracies. The outliers above the respective theoreti- 
cal lines for each individual pulsar are described below: 

1 . Radio ephemerides extrapolated over long time intervals ap- 
pear to be unreliable. Therefore, we exclude the following 
observations in the final calculation regarding the relative 
timing accuracy: Vela pulsar: 015395140; PSR B 1509-58: 
01281 2040 1 ; PSR B0540-69: 04 1 3 1 8020 1 , 04 1 3 1 8030 1 . 

2. The 6P approximation given in Eq. |4]was found to be un- 
reliable in some cases (e.g. for PSR B 1509-58, observation: 
0312590101 and the Vela pulsar, observations: 01 1 1080101, 
0111080201). As the Vela pulsar is quite active, we may 
have under-estimated the error by simply extrapolating 
the ephemeris, however, the same can not be said about 
PSR B 1509-58 which is one of the most stable young pul- 
sars known. 

As seen in Fig. [6] there is no obvious change in the relative 
timing accuracy of the EPIC-pn camera over its lifetime. 

The results for the Crab pulsar alone are shown in Fig. [7] As 
expected, there is a tendency towards smaller uncertainties for 
longer observations. 

For a quantitative measure of the timing accuracy the stan- 
dard deviation for the AP/P distribution (shown in Fig. [6]l was 



10 -4 


I I I 1 

+ PSRJ0537- 


69 










1 1 1 




i 


A Crab 

* PSRB0540 

ovelo 

□ PSRB1509 


-58 
















i rv-5 
10 


X PSRB1055 

□ 


-52 


























O 








10' 6 


o 










□ 






- 


lO" 7 


- H 












A 




A 


lO" 8 


~A 


A 
A 


A 




A 
\* 


A 


A 


A 


A 






A 


A 








& 


A 




lO" 9 




A 








A 


A A 


A 


10 10 


A 

, , , i 


+ 


, , r .1 


A 






A 




A 



1000 2000 3000 

Epoch [MJD-51000] 



Fig. 6. The relative timing accuracy of XMM-Newton 's EPIC-pn 
camera: |AP/P| for all pulsars as a function of date [MJD]. 

used. Fitting the distributions with a Gaussian normal distribu- 
tion, we found a standard deviation of 7 x 10 9 for all the pulsars 
together (including the Crab pulsar) and 5 x 10~ 9 for the Crab 
pulsar alone. While the distribution for Crab pulsar is centred at 
zero (within uncertainties) the mean value of the distribution for 
all the pulsars combined is slightly offset, in the sense that the 
X-ray period is slightly shorter on average than the radio period. 
Thus, the relative deviation of the observed pulse period with 
respect to the most accurate radio data available is AP/P < 10~ 8 . 

5. Absolute timing accuracy of XMM-Newton 

The Crab pulsar shows a shift of -306 jus (shown in Fig. [8]) be- 
tween the peak of the first X-ray pulse with respect to the radio 
peak. We here by confirm the simil ar results of o ther missions 
like C handra dTennant et all 1200 ll) INTEGRAL dKuiper et all 
120031) . and RXTE dRots et all |2004|) . The XMM-Newton values 
("stars" in Fig. [5J show a considerable scatter with a standard de- 
viation of 48 /is. The formal error on the mean value of -306 fis 
is ±9fis. 

The scatter found is consistent with the previously deter- 
mined maxim um integrated error for the time correlation of less 
than 100 //s dKirsch et all l2005h . The original requirement for 
an absolute timing accuracy of 1 ms for XMM-Newton, defined 
before launch, has clearly been reached and even improved on 
by at least a factor of 20. 

This scatter is likely to be due to uncertainties in the time 
correlation process since the phase of the main peak can be mea- 
sured with an ac curacy of yus. Upper limits for these processes 
were reported by Kirsc h et aD (12004 who gave a detailed de- 
scription of all kinds of instrumental delays considered while 
converting between observing time and UTC time and estimated 
the spacecraft clock error to be ~ 1 1 ps, the uncertainty in ground 
station delays to be ~ 5 fis, the interpolation errors to be ~ 10/ks, 
the error between latching observing time and the start of frame 
transmission as ~9/ks, and the uncertainties in the spacecraft 
orbit ephemeris to be ~30/is. All these errors will be random 
for our data, and hence the fluctuations observed. The 48 [is 1 cr 
scatter measured with respect to the mean may then be attributed 
uniquely to the above errors and no other systematic effect. This 

6 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0045-l-0.pdf 
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Fig. 7. Relative timing using Crab monitoring data: upper panel 
as a function of observing time, lower panel as a function of 
MJD. These plots are regularly updated i n the EPIC Calibration 
Status Document dGuainazzi et all 1201 lb using the routine cali- 
bration observations of the Crab. 



value can then be considered to be the minimal significant time 
separation between two arrival times to be considered indepen- 
dent. 

From the initial 38 Crab observations, 32 were considered 
for the absolute analysis. Six of the Crab observations were ex- 
cluded for the reasons given below: 

- Observation 0122330801: was early in the mission (rev. 56) 
and appears to have problems in the time correlation that can 
no longer be recovered. 

- Observation 0160960201: too much data had to be excluded 
due to time jumps which caused a dramatic reduction in the 
number of counts (rev 698). 

- Observations 0160960401 and 0160960601: these corre- 
spond to rev. 874 which shows the X-ray pulse peak dis- 
placed from the expected radio position. This is likely to 
have been caused by a glitch shortly before the XMM- 
Newton observations. This offset is more dramatic in the sec- 
ond observation, which has poorer statistics due to a shorter 
exposure time. None of these observations fall into the range 
shown in Fig. [8] 



- Observations 0412590601 and 0412590701: these corre- 
spond to rev. 1325. The reason for the offsets is unclear. They 
may be due to a small, non-reported glitch or an anomalously 
large ground segment error, but because of the uncertainty, 
we excluded these observations when determining the abso- 
lute timing precision. 



It was found that some observations presented pulse profiles 
with an excess of counts in the interpulse region of the Crab pro- 
file. Numerical simulations have been used to study the effect 
that this excess could have caused in determining the peak of the 
X-ray profile and thus, in determining the difference in phase 
between the X-ray and the radio. Using the typical 0.2-15.0 keV 
Crab profile as the input, 10000 light curves were created us- 
ing Monte Carlo simulations. The strength of both peaks in the 
pulse profile (keeping the ratio between them constant) as well 
as the strength of the interpulse were selected as input parame- 
ters. As shown in Fig. IC.ll in Appendix C, the Moffat fit, used 
to determine the phase of the main peak, can be used to reliably 
fit different tails. To fit the excess in the interpulse we added 
a Lorentzian function to the pulse profile in order to take into 
account the excess in the tail of the main peak. The secondary 
peak was omitted in the fit. No phase shift was found in any of 
the models tried, which implies that the strength of the interpulse 
region plays no role in determining the phase of the peak of the 
profile and thus all 32 of the retained Crab observations could be 
used to derive the absolute timing accuracy reliably. 
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Fig. 8. Absolute timing as a function of MJD using Crab-pulsar 
monitoring data. Shown are the offsets of the main X-ray peak 
with respect to the main radio peak in the Crab pulse profile in 
time units (left scale) and phase units (right scale). The colored 
lines give the mean values for XMM-Newton (solid black line, 
this work), Chandra (green dashed line), RXTE (red dot line) and 
INTEGRAL (blue dashed-dot line) respectively, all taken from 
the literature. The yellow area indicates the standard deviation of 
the XMM-Newton data points. The superscript numbers near to 
each XMM-Newton data point give the XMM-Newton revolution 
in which the observation was carried out. 
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6. Discussion and Conclusion 

The Crab pulsar has been used by many missions as a 
calibration source for timing accuracy (iKuiper et al.l (|2003 | ) . 
[ Rots et all (l2004.IOosterbroek et all (l2008l) . lAbdo et al.l (120 id) . 
iMolkov et al. Il2010h . The XMM-Newton observatory began ob- 
serving the Crab pulsar during its earliest orbits, monitoring its 
X-ray pulsation with high time resolution. 38 Crab observations 
spread over 10 years have been analysed in this paper (from rev- 
olution 56 until revolution 1687). Measurements of the period 
were made with an accuracy of ~10~ n s. A relative timing ac- 
curacy smaller than 10~ 8 and stable with time was established 
for the EPIC-pn camera. This result was achieved by comparing 
our X-ray measurements of the Crab pulsar with high precision 
radio measurements at each corresponding epoch. Five isolated 
pulsars showing a wide range of periods and completely differ- 
ent pulse profiles (PSR J0537-69, PSR B0540-69, Vela pulsar, 
PSR B 1509-58 and PSR B 1055-52) were analysed to comple- 
ment the study of the relative timing accuracy, confirming the 
results obtained with the Crab pulsar. 

For the case of PSR B0540-69 a long term phase-coheren t 
study of its period was reported by Livingston e et al.l d2005b). 
Due to its stability we considered it a good candidate to use for 
an extrapolation over a long time period. As shown in Table [2] 
the long extrapolation made in two of the three observations of 
PSR B0540-69 show poor results suggesting that a small glitch 
between the ephemeris and our observation may have occurred, 
rendering this pulsar less stable than anticipated. 

An improved algorithm to detect and correct sporadic 
"jumps" in the flow of t he photon arrival tim es has been imple- 
mented with SASv8.0 (Guai nazzi et all 1201 1 R This method is 
based on a more accurate determination of the frame times for all 
pn modes and on a correction of frame time drifts due to temper- 
ature variation and aging of the on-board clock (Frevberg et al., 
120051) . The total reduction of the rate per 100 ks of observation 
affected by residual uncorrected time jumps for all pn instrumen- 
tal modes dropped from 2.8 before the improved algorithm to 0.3 
once it was implemented. 

For the absolute timing analysis, only Crab pulsar observa- 
tions have been analysed since a high number of stable obser- 
vations need to be considered to provide a reliable result. We 
have considered the phase of the first (main) peak of the X-ray 
profile and measured the phase difference with respect to the 
corresponding peak of the radio profile. Considering 32 of 38 
Crab EPIC-pn observations (0.2-12 keV energy range) analysed 
in this paper, we confirmed previous results demonstrating that 
the first X-ray peak from the Crab pulsar leads the radio peak by 
306 + 9/is (statistical error) with +48 /vs (1 cr) scatter. This error 
is similar to the Ground Segment accuracy and defines the abso- 
lute timing accuracy of the instrument. T he observed shift is con- 
sistent within 1 cr with th ose presented bylKuiper et al.l d2003l) us- 
ing INTEGRAL and bv lRots et al] d2004h using RXTE, as shown 
in Fig. E 

A systematic comparison of our measurements in the X-ray 
band with respect to other accurate measurements carried out 
in different energy bands from earlier observations in the opti- 
cal, X-ray, and y-ray parts of the spectrum are shown in Fig. 
Differences in the shifts observed over 7 decades in energy are 
marginal with an average value of 284.4/ys. It is important to 
note that the large error bars quoted in the X-ray band for XMM- 
Newton and RXTE include systematic errors from the radio mea- 
surements, carried out at the Jodrell Bank Observatory. 
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Fig. 9. The peak pulse lead time (//s) of various observations are 
plotted against energy (keV) in optical, X-rays and y-ray energy 
bands and a constant model is fitted to the data, which is found 
to be 284.4 fis. The data points shown above cover over 7 orders 
of magnitude in energy and come from d ifferent obse r vation s 
and experi ments mentioned on the graph (IKuiper et ail (12003 1 ) . 
[ Rots et al J (120041). lOosterbroek et al.l (120081) . lAbdo et al.l (120 1 01) , 
Mol kov et al. I d2010h ). Only the central energy corresponding to 
an individual observation was plotted along with corresponding 
lead time measured within the observed energy band. 



The origin of the electromagnetic radiation emitted from pul- 
sars is still unclear. Several models have been proposed to ex- 
plain the origin of the high energy radiation based on differ- 
ent regions of acceleration in the pulsar magnetosphere, such 
as the polar cap, the slot gap and the outer gap models (e.g . 
Harding et al., 1978; Arons & Scharlemanh , 1 19791: ICheng et all 
1986; Zhang & Harding! 120001; lHarding & MuslimovL \20lW . 
The radio emission model is an empirical one and the radiation 
is usually assumed to come from a core beam centered on the 
magnetic axis and one or more hollow cones surrounding the 
core (e.g. Rankin, 1983). The estimated average delay between 
the emission from differing wavelengths is therefore very sig- 
nificant. It implies that the site of radio production is distinctly 
different from that of the non-radio emission. The difference in 
phase between the radio and the X-ray radiation is about 0.008, 
or three degrees in phase angle. This time delay of about 300 fis 
most naturally implies that emission regions differ in position by 
about 90 km between radio and X-rays energy bands in a sim- 
plistic geometrical model neglecting any relativistic effects, with 
the radio being emitted from closer to the surface of the neutron 
star. Such high time resolution, high precision absolute timing, 
multiwavelength observations are therefore essential for under- 
standing the origin of the pulsar emission. 
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Appendix A: XMM-Newton and its fast EPIC-pn 
timing modes 

EPIC is capable of providing moderate energy resolution spec- 
troscopy in the energy band from 0.2 to 15 keV for as many as 
several hundred sources in its 30' field-of-view. The EPIC cam- 
eras can be operated in different observational modes related to 
different readout procedures. Detailed descriptions of the vari- 
ou s readout modes of EPIC-pn and their limita tion s are given 
bvlKendziorraet al.1 (1 19991) . iKuster et a l. ( 1999) and Ness et al.1 
(120 KM 

The EPIC-pn camera provides the highest time resolution 
in its fast Timing and Burst modes (Timing mode: 29.52yus, 
Burst mode: 7 yus) and moderate energy resolution (E/dE = 10- 
50) in t he 0.2-15 k e V ene rgy band. The pile-up limit (see Sec. 
3.3.9 of iNess et al.l (12010b ') for a point source is 800 counts s _1 
(85mCrab) in Timing mode and 60000 counts s _1 (6.3 Crab) in 
Burst mode. Thus, the observations of the Crab suffer from pile- 
up only in Timing Mode, such that spectral analysis of the Crab 
can only be carried out accurately in Burst mode. However for 
timing purposes the effect of pile up can be neglected in the 
Timing mode. 




duced in CCD coordinates using RAWX and RAWY, which are 
simply the pixel co-ordinates, where each pn pixel is 4.1x4.1" 
aside. The source appears as a stripe in the CCD RAWY direc- 
tion. Source extraction regions in both modes will therefore al- 
ways be boxes. 

Operated in Timing mode EPIC-pn data show a bright line 
in the RAWX direction at RAWY = 19, that is related to a fea- 
ture in the on-board clock sequence. In the clocking scheme of 
the Timing mode 10 lines are shifted to the CAMEXs (CMOS 
Amplierand Multiplexer Chip) and then read out as one so called 
macro line, such that the integration time for a normal macro line 
is 29.52yus. Within a frame time, 200 macro lines are read (cor- 
responding to 2000 physical CCD lines). During the first CCD 
readout the first macro line contains only one CCD line and is 
set to bad. However, the integration time during the readout of 
the second macro-line is 29.52 + 23 //s due to electronic imple- 
mentation of the sequencer. Therefore the integration time for a 
point source at RAWY = 189 is a factor 1.8 higher than for all 
other macro line and macro line 19 receives a factor 1.8 higher 
flux from the point source. The feature only shows up for bright 
point sources. 

There is no effect on the scientific quality of the data as long 
as the integration time for spectra and light curve s is higher than 
the fr ame time in Timing mode (5.96464ms; dFrevberg et al.l 
l2005lPl Caution should be used for pulse phase spectroscopy 
with bin sizes below the frame time (5.96464ms), but only if 
the pulse period is a multiple of the frame time. 
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Fig. A.l. Upper: RAWX and RAWY image of a Crab observation 
in Burst mode. The rows 180-200 are not transmitted. Bottom: 
RAWX and RAWY image of a Crab observation in Timing mode. 



As shown in Fig. IA.ll in both timing and burst modes, EPIC- 
pn loses spatial resolution in the shift-direction. In Timing mode, 
this is because 10 lines of events are fast-shifted towards the an- 
odes and then the integrated signal is read out as one line. In 
Burst Mode, it is because 200 lines are fast-shifted within 14.4 fis 
while still accumulating information from the source. The stored 
information is then read out as normal, where the last 20 lines 
have to be deleted due to contamination by the source during 
the readout. The CCD is then erased by a fast shift of 200 lines, 
and immediately after that the next Burst readout cycle starts. 
Moreover, the lifetime in Burst mode is only 3% and therefore, 
the use of this mode has been limited to observations of very 
bright sources such as the Crab or X-ray transients. For our anal- 
ysis we use mainly Timing and Burst mode observations. The 
images seen in Fig. lA.ll for the Timing and Burst mode are pro- 

8 http://xmm.esac.esa.int/external/xmm_user_support/ 

documentation/uhb/XMM_UHB . pdf 9 http : //xmm2 . esac . esa . int/docs/documents/CAL-TN-088 1 . pdf 
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Appendix B: Treatment of uncertainties and 
reliability of radio extrapolations 

The x 2 distribution obtained from the period search can be ap- 
proximated by a triangle where the maximum corresponds to the 
true period Po and the points Pi and P2 where the legs of the tri- 
angle meet the level of constant^ 2 defining the total width of the 
X 2 distribution. For a pulse profile with a small single peak, Pi 
and P2 can be calculated using Eq. IB.ll where T b s is the elapsed 
observational time and N per is the number of pulse periods in this 
time. 



I obs 



N peI + 1 



Pi 



I obs 



1 



where N D 



T'obs 



(B.l) 



For a triangular function the Full Width Half Maximum 
(FWHM) is equal to (P2 - P\)/2 and can be expressed as in 
Eq. IB. 21 as a function of the period and the elapsed observation 
time. 

Pi - P\ P 2 
-=— — - => FWHM = 

2 r bs 
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Fig. B.l. Sample^ 2 distributions for one observation of each of 
the studied pulsars. 

An expected FWHM of the x 2 distribution can be estimated 
using Eq. IB. 21 A comparison of such estimations (lines) and 
the measured FWHMs from all the observations is shown in 
Fig. IB. 21 All values were normalised using the pulsar period to 
be able to present all the pulsars on the same diagram. All mea- 
sured values of FWHM/P are about a factor of 3 smaller than 
those predicted for the Crab and Vela pulsars. In the other four 
pulsars the ratio between the measured and predicted values is 
~ 1.3. This would suggest that this approximation works better 
in single peaked pulse profiles. 

Empirically, two periods can be considered completely inde- 
pendent from each other when their difference is at least P 2 / T 
(one Independent Fourier Space, IFS (Ide Ona WilhelmiL l2003l) 1. 
which is identical to the FWHM definition in Eq. IB. 21 One IFS 
can be seen, then, as the delta in pulse period which will smear 
a perfect pulse profile to a flat profile when folded over the com- 
plete observing time r b S . This approach is quite conservative 
and smaller changes than one IFS in period can be easily seen. 
We have found that a rough estimate of the uncertainty in the 




10* 

Time Base [s] 



Fig. B.2. Comparison of the predicted FWHM of the x" distri- 
butions (lines) and the observed ones (symbols). All values were 
normalised using the pulsar period to be able to present all the 
pulsars on the same diagram. 



measured period can be found by dividing the FWHM by the 
number of phase bins used to construct the pulse profile (degrees 
of freedom, see Table [3}. Thus, two periods will be considered 
different when the pulse profile is smeared by one bin instead 
of one whole phase. The error on the X-ray period can then be 
written as shown in Eq. IB. 31 



5P 



FWHM 

dof 



(B.3) 



Besides providing a good estimate of the error on the period, 
the Independent Fourier Space approach can also provide a good 
indication of how reliable the extrapolation (or interpolation in 
the case of the Crab pulsar) of a period can be. Since we defined 
the relative timing accuracy in Sect. I3.1l based on the reference 
period (normally obtained from radio observations) at the time 
of the XMM-Newton observation, it is critical to understand how 
reliable, and in some form, how accurate, this parameter really 
is. For clarification, and due to the huge amount of data available 
we will focus on the Crab pulsar only. However the same applies 
to all the pulsars studied in this paper. 

Eq. lB.3l establishes that two periods are completely different 
if the pulse profile is smeared by one bin. By studying the phase 
smear, we can determine whether the pulse profile has been af- 
fected by a glitch or whether extrapolating the ephemeris over 
(long) time periods leads to inaccuracies. 

If a simple period evolution with time (including the sec- 
ond derivative) is assumed, the phase smear is then defined as in 
Eq.lBH 



Phase Smear 



(P 



extrap 



Po)T, 



obs 



^extrap X Pq 



(B.4) 



where P ex trap is the extrapolated period at the time To, Po is the 
actual period at that time and T bs is the exposure time of the 
observation. 

Using radio data from the Jodrell Bank Observatory, the 
phase smear versus time of extrapolated periods is shown in 
Fig. IB. 31 For the relative timing analysis 100 phase bins have 
been used and therefore a limit of 1 % smearing is imposed by 
the criteria described above. Extrapolating the period, the Crab 
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pulsar would reach that limit in 4 months. Considering that the 
Jodrell Bank Observatory provides an updated ephemeris every 
month, the actual smearing effect will be much lower (~ 0.1%) 
and other properties such as timing noise will not affect our rela- 
tive timing analysis. For the Crab, we actually used interpolation 
rather than extrapolation, see Sect. 13. II so the phase smear was 
further minimised (~ 0.09%). 
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Fig. B.3. Phase smear of the Crab pulsar versus time. The tim- 
ing noise seems to be a prominent feature for the Crab pulsar 
and therefore extrapolation over periods of 4 months will already 
reach our 1% limit (dashed line). 
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Appendix C: The Moffat function 

The Moffat functio n is a mo dified Lorentzian with a variable 
power law index (Moffal fl969h . In Fig. O the behaviour of the 
function is shown as a function of its parameters. The function 
presents different tails on each side of the maximum which fit the 
main pulse of the Crab profile better than a normal Lorentzian or 
Gaussian function. The explicit formula of the Moffat function 
is the following: 



y = 777 — 7-777-7 — rrr +A 4 +A 5^ (C.i) 

(((x-Ai)/A 2 ) 2 + I) 71 ' 
The different parameters represent: 

Aq: normalization 
A\\ Peak Centroid 
A 2 : HWHM 
A 3 : Moffat index 
A4: offset 
A5: slope 

The variation in the shape of the Moffat function for different 
values of the important parameters is shown in Fig. IC.ll Upper 
left: A2 changes from 0.02 to 0.06; upper right: A 5 changes from 
500 to 1400; lower left: A\ changes from 0.4 to 1.12; and lower 
right: A3 changes from 1.0 to 2.80. 
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